Control system for extraction of hydrocarbons from underground deposits

ABSTRACT

A control system for use in extracting hydrocarbons from an underground deposit is disclosed that comprises an electromagnetic heating system and a processor. The electromagnetic heating system is configured to heat the underground deposit to facilitate fluid flow of a resource for extraction from the underground deposit. The processor is configured to control the electromagnetic heating system in response to temperature data and pressure data for the underground deposit. The processor correlates the temperature data and pressure data with predetermined water phase characteristics to control the electromagnetic heating system so that substantially all water in the underground deposit is maintained in a liquid state. The control system may also generate voxel data corresponding to spatial characteristics of the underground deposit. The spatial characteristics may be presented as a map on a display.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

[Not Applicable]

BACKGROUND OF THE INVENTION

As the world's standard crude oil reserves are exhausted, and the continued demand for oil causes oil prices to rise, oil producers are trying to increase their output of hydrocarbons from bituminous ore, oil sands, tar sands, and heavy oil deposits. These materials are often found in naturally occurring mixtures of sand or clay. Because of the high viscosity of bituminous ore, oil sands, oil shale, tar sands, and heavy oil, the drilling and refinement methods used in extracting standard crude oil are typically not available. Therefore, recovery of oil from these deposits requires heating to separate hydrocarbons from other geological materials and to preserve hydrocarbons at temperatures at which they will flow. Steam is typically used to provide this heat, although electric and radio frequency heating is sometimes employed. The heating and processing can take place in-situ, or in another location after strip mining the deposits.

During in-situ processing, it is difficult to sense the real-time conditions in the deposit and the well bore. Uncertainty often exists about whether the hydrocarbons are indeed flowing. Valuable time and heat energy are wasted in unsuccessful attempts at in-situ processing when low formation permeability conditions prevent diffusion of steam and limit heating, or when steam and heat move away from a zone that is targeted for heating through formation fractures or through high permeability materials.

“Well logging” may be used to check and record subsurface conditions. Such well logging involves an examination of core samples, and moving sensors up or down the well bores. Sensors are used to measure, for example, electrical resistance, acoustic properties, natural radioactivity, density of the formation surrounding the well, and the like. However, these measurements do not produce a near real-time picture of conditions in the formation. They produce only a static and partial picture of such conditions.

SUMMARY OF THE INVENTION

A control system for use in extracting hydrocarbons from an underground deposit is disclosed that comprises an electromagnetic heating system and a processor. The electromagnetic heating system is configured to heat the underground deposit to facilitate fluid flow of a resource for extraction from the underground deposit. The processor is configured to control the electromagnetic heating system in response to temperature data and pressure data for the underground deposit. The processor correlates the temperature data and pressure data with predetermined water phase characteristics to control the electromagnetic heating system so that substantially all water in the underground deposit is maintained in a liquid state. The control system may also generate voxel data corresponding to spatial characteristics of the underground deposit. The spatial characteristics may be presented to a user as a map on a display.

Other aspects of the invention will be apparent from this disclosure.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows one embodiment of a steam assisted gravity drainage system for in-situ extraction of hydrocarbons from an underground resource, where the extraction system includes a control system.

FIG. 2 is a block diagram of an embodiment of the system shown in FIG. 1.

FIG. 3 is a block diagram of elements and operations included in and executed by the system shown in FIG. 2.

FIGS. 4 and 5 are Mollier charts showing temperature on the x-axis and pressure on the y-axis for different water phases.

FIG. 6 shows one embodiment of a kill switch.

FIG. 7 illustrates the impedance encountered at various portions of the embodiment of the system shown in FIG. 2.

FIG. 8 shows one embodiment of a display system 800 that may be used to display a map of the characteristics of the underground deposit.

FIG. 9 is one embodiment of a voxel manipulation system 900 that may be used in the system shown in FIG. 8.

FIG. 10 is a table of coefficients that may be used to generate dielectric voxels at various operating frequencies of the electromagnetic heating system.

FIGS. 11 and 12 are examples of maps that may be displayed using the system shown in FIG. 8.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 illustrates one embodiment of a steam assisted gravity drainage system 15 for in-situ extraction of hydrocarbons from an underground deposit 20. As shown, steam assisted gravity drainage system 15 includes a control system 10 that controls the operation of the various components of the system. Steam and other suitable susceptors for the steam assisted gravity drainage system 15 are provided by a steam and susceptor system 25. The steam and susceptors are injected from the steam and susceptor system 25 through a delivery pipe 35 and into the underground deposit 20. The steam and susceptor heat the materials in the underground deposit 20 to a temperature at which hydrocarbons in the formation, condensed steam, and possibly other liquids drain down to a recovery pipe 40 where they are pumped out as recovered resources into storage tanks 42 for storage and processing. Although control system 10 is described in the context of the steam assisted gravity drainage system 15, it may be employed in any type of hydrocarbon extraction system.

Electromagnetic energy may be used along with the heated mixture of steam and susceptor to augment the heating of the underground deposit 20. To this end, steam assisted gravity drainage system 15 includes an electromagnetic heating system 45 that provides electromagnetic energy to underground deposit 20. Electromagnetic energy is provided from an RF source 47 to an antenna 50. Electric and magnetic fields are generated at antenna 50 to heat the materials in underground deposit 20. Magnetic fields generated at antenna 50 heat materials in underground deposit 20 that exhibit magnetic dissipation while conductive materials are heated by eddy currents that can be induced by both magnetic and electric fields. Such materials are efficiently heated by RF energy when the strongest fields generated at antenna 50 are matched with the materials in underground deposit 20. For example, conductive materials such as water and, particularly, water mixed with sodium hydroxide, are efficiently heated by eddy current produced by a magnetic field at the antenna 50. Materials that are not electrically conductive but that exhibit dielectric dissipation are effectively heated at antenna 50 by RF electric fields.

Mixtures of materials (e.g., steam and susceptors) that are susceptible to electromagnetic heating may be injected through delivery pipe 35 and into underground deposit 20 from the surface. The resulting heated mixture, including the hydrocarbons, are all electromagnetic the heated to a point at which they are liquefied and pumped from underground deposit 20 through recovery pipe 40. Copending applications assigned to Harris having docket numbers 20478US01 and 20483US01 disclose heating of hydrocarbons by mixing hydrocarbons with materials that are strongly susceptible to electromagnetic heating. Both applications are included herein by reference.

Various facets of the steam assisted gravity drainage system 15 may be controlled by the multiphysics control system 55 using sensed/calculated data 60. The sensed/calculated data 60 may include, for example, real-time temperature and pressure data 62, resource phase state data 64, VSWR values 66, recovered resource temperature data 68, seismic acoustic data 70, dead reckoning data 72, and recovered resource chemistry data 74.

The sensed/calculated data 60 represents various characteristics of the underground deposit 20. More particularly, the resource phase state data 64 may be used to determine whether the water and/or other resources in the underground deposit 20 are in gas or liquid states. The VSWR values 66 provide an indication of the extent to which the output impedance of electromagnetic heating system 45 matches the input impedance of the antenna 50 and the underground deposit 20. Minimizing the VSWR values 66 ensures proper power transfer from the electromagnetic heating system 45 to the underground deposit 20. The VSWR values 66 may be determined from pulse reflectometry. The recovered resource temperature data 68 correspond to the temperature of the resources in the underground deposit 20. The seismic acoustic data 70 correspond to data obtained from, for example, echograms of the underground deposit 20. The dead reckoning data 72 corresponds to such things as water injected to date, watts used to date, and similar consumption of materials and energy used for the extraction of the resource from the underground deposit 20. The recovered resource chemistry data 74 may include data such as the molecular weight, water-to-oil ratio, recovered resource water phase, and other characteristics of the underground deposit 20.

Among other things, the multiphysics control system 55 uses the sensed/calculated data 60 to control provision of electromagnetic energy, steam, and susceptor to the underground deposit 20. In the system shown in FIG. 1, multiphysics control system 55 controls the temperature and composition of the steam and susceptor components of susceptor system 25 as well as the amount of the steam and susceptor provided through delivery pipe 35 to underground deposit 20. The steam and susceptor provided through delivery pipe 35 may be controlled, at least in part, through manipulation of the elements of the motor and valve system 80. The multiphysics control system 55 may also use the sensed/calculated data 60 to control parameters of the RF energy provided from electromagnetic heating system 45 to the underground deposit 20. As will be explained in further detail below, multiphysics control system 55 controls the parameters of the applied RF energy based, at least in part, on the phase of the water in the underground deposit 20. The water phase may be determined using one or more of the sensed/calculated data 60.

FIG. 2 is a diagram of a further embodiment of steam assisted gravity drainage system 15. In this embodiment, the steam assisted gravity drainage system 15 employs a multiphysics control system 55 that receives data representing the characteristics of the underground deposit 20 and uses that data to control operation of the various components of the system 15. Some of this data may be obtained from an analysis of extracted resources stored in storage tanks 42. Such recovered resources are extracted through recovery pipe 40 and passed through a gas and liquid separator 130 before being provided to the storage tanks 42. The storage tanks 42 may include an oil storage tank 115, a gas storage tank 120, and a further storage tank 125, where storage tank 125 may be used to store miscellaneous constituents of the recovered resource. The extracted resource analysis may be executed, for example, as an off-line process.

The data used by the multiphysics control system 55 is not limited to data obtained by the extracted resource analysis. Rather, data may also be obtained on a substantially real-time basis from sensors placed at strategic locations in steam assisted gravity drainage system 15 (see in-situ sensor network 310 of FIG. 3). Such data may include, for example, real-time temperature and pressure measurements of the underground deposit 20.

The characteristics of underground deposit 20 provided to multiphysics control system 55 in FIG. 2 include the conductivity 135 of the resources in underground deposit 20, the temperature 140 of the resources and underground deposit 20, the permeability 145 of the resources in of underground deposit 20, and the pressure 150 in underground deposit 20. The multiphysics control system 55 uses this data to exercise control over the elements of steam assisted gravity drainage system 15 that apply electromagnetic energy, steam, and susceptor to underground deposit 20. One of the goals of such control is to increase the efficiency of the energy transfer from the system 15 to the underground deposit 20. This may be accomplished, for example, by controlling the parameters of the RF power applied to underground deposit 20 as well as controlling the composition and injection of the steam and susceptor to underground deposit 20.

As shown in FIG. 2, the system 15 provides steam and susceptor to the underground deposit 20 using a fluid supply system 155. The fluid supply system 155 may employ storage tanks 160 including, for example, a sodium chloride supply tank 161, a water supply tank 162, and a caustic soda supply tank 163. The multiphysics control system 55 manipulates the composition of the mixture provided to underground deposit 20 by regulating the flow of liquid provided from the storage tanks 160 to a susceptor mixing and injection system 170. In this example, multiphysics control system 55 may regulate this flow by providing the appropriate control signals to a valve controller 165 disposed between the storage tanks 160 and the susceptor mixing and injection system 170. The multiphysics control system 55 operates the valve controller 165 based on conductivity 135, temperature 140, permeability 145, and/or pressure 150 to provide a mixture having the desired composition. Once the desired mixture is obtained, it is provided to the susceptor mixing and injection system 170 and, therefrom to delivery pipe 35 through fluid pump 175.

Electromagnetic energy is provided to the underground deposit 20 using the electromagnetic heating system 45 to heat the resources in the underground deposit 20. In the embodiment shown in FIG. 2, the electromagnetic heating system 45 includes an RF signal generator 185 having its frequency, power, and waveform modulation controlled by the multiphysics control system 55. The RF signal from electromagnetic heating system 45 is coupled to the input of an impedance matching network 190, which is used by the multiphysics control system 55 to calculate the VSWR values 66. The multiphysics control system 55 controls the impedance matching network 190 to minimize the VSWR values 66 thereby maximizing power transfer between the RF signal generator 185, antenna 50, and underground deposit 20. As shown in FIG. 2, the output of the impedance matching network 190 may be provided to the input of a directional coupler 195 having a first terminal coupled to a load 200 and a second terminal coupled to the antenna 50. The electromagnetic heating system 45 may also include RF chokes 205 that are positioned to prevent undesired electromagnetic energy leakage from the electromagnetic heating system 45 to other portions of the system 15.

The multiphysics control system 55 is configured to control electromagnetic heating system 45 in response to the temperature and pressure of the underground deposit 20. Referring to FIG. 1, the temperature and pressure data 62 may include real-time temperature and pressure measurements obtained from sensors positioned at strategic locations in system 15 and/or temperature and pressure data 62 derived through the extracted resource analysis.

The sensed temperature and pressure of the well 340 is correlated by the multiphysics control system 55 with predetermined water phase characteristics to control the operation of the electromagnetic heating system 45 so that substantially all water in the underground deposit 20 is maintained in a liquid state. N one example, multiphysics control system 55 is configured to use sensed temperature measurements and sensed pressure measurements in conjunction with Mollier water phase temperature and pressure data to determine the phase of the water in underground deposit 20. The multiphysics control system 55 may use this water phase determination to provide signals to electromagnetic heating system 45 that control waveform parameters such as frequency, modulation, and power. Further, if attempts to maintain the water in a liquid phase using the waveform parameters are unsuccessful, electromagnetic heating system 45 may be disconnected so that it no longer provides any electromagnetic energy to the underground deposit 20.

The multiphysics control system 55 may further correlate the temperature and pressure characteristics of underground deposit 20 with its permittivity and conductivity characteristics. For example, the temperature and pressure measurements may be correlated with Mollier data, which, in turn, may be used to determine the permittivity and conductivity of resources in underground deposit 20. The multiphysics control system 55 may use the permittivity and conductivity measurements obtained in this manner as alternative or additional data that, in turn, may be used to control the electrical parameters of the electromagnetic heating system 45.

FIG. 3 is a block diagram of elements and operations included in and executed by the system shown in FIG. 2. In this example, various components of system 15 provide input data to multiphysics control system 55, which, in turn, is used to generate output control signals. Among the output control signals are those that are used to apply electromagnetic energy, steam, and susceptor to the underground deposit 20.

The multiphysics control system 55 may begin its control operations using an a priori model of system 15 and underground deposit 20. The a priori model may be generated using well logging information obtained before the well 340 is drilled. Such well logging information may include electrical resistance maps, alpha particle backscatter information, and the like.

At operation 305, multiphysics control system 55 derives well characteristics using the well behavior model. The initial a priori well behavior model may include values for the RF power and frequency of electromagnetic heating system 45, steam pressure, and resource spatial distribution. These values may be derived from well logging and other pre-drilling information. Further, the a priori model may be used as a means to provide initial temperature and pressure values for underground deposit 20. Still further, the a priori model may be used as a means to provide initial bulk permittivity (∈) and bulk conductivity (σ) of underground deposit 20.

As the well 340 is used over time, new well data becomes available to multiphysics control system 55. This new well data, in turn, is used to update the a priori model to a more accurate working well behavior model. The updated data may include resource properties obtained during a recovered resource analysis. Additionally, or in the alternative, the well behavior model may be updated using real-time sensed well parameters obtained from an in-situ sensor network 310. The well behavior model may be updated at predetermined intervals or whenever an updated version of the well behavior model is desired.

The frequency, power, modulation, and impedance parameters are used by multiphysics control system 55 at operation 307 to determine whether there is an adequate impedance match between the electromagnetic heating system 45 and underground deposit 20. Time domain reflectrometry may be used in this determination. To this end, the impedance matching network 190 may include a time domain reflectometer that is connected to delivery pipe 35 and recovery pipe 40 through directional coupler 195. Delivery pipe 35 and recovery pipe 40 may serve as antennas for the time domain reflectometer and form a transmission line of the “open wire” type. During the extraction process, the RF signal generator 185 may transmit a short rise time RF pulse to delivery pipe 35. The characteristics of the reflected pulse are a function of the impedance encountered by the pulse and may be used to determine the VSWR values.

The multiphysics control system 55 performs different operations depending on whether the impedance match is satisfactory. If it is satisfactory, the multiphysics control system 55 controls RF signal generator 185 at operation 315 to provide electromagnetic power through kill switch 320 for ultimate provision to underground deposit 20. Some of the electrical parameters that may be used at operation 315 to control RF signal generator 185 include waveform frequency (f), power (P_(RF)), modulation frequency and type, and, as noted above, impedance (Z). However, when the impedance match is inadequate, multiphysics control system 55 adjusts the impedance of the impedance matching network 190 at operation 325 to minimize the VSWR thereby maximizing power transfer to the underground deposit 20 and.

The multiphysics control system 55 decides at operation 330 whether steam has formed in underground deposit 20. To this end, multiphysics control system 55 may receive real-time temperature (t) and pressure (P) data from in-situ sensor network 310 and correlate this data with, for example, Mollier water phase data. If steam is detected, multiphysics control system 55 may adjust the parameters of electromagnetic heating system 45 at operation 315 to prevent further steam generation in underground deposit 20. Alternatively, multiphysics control system 55 may operate kill switch 320 to effectively disconnect electromagnetic heating system 45 from underground deposit 20 to prevent further application of electromagnetic energy.

If steam is not detected, multiphysics control system 55 controls the boiler at operation 345. Steam from the boiler is provided to a boiler steam supply at operation 350. The steam and solvent are mixed at operation 355 before being injected into underground deposit 20 at operation 360.

The multiphysics control system 55 may also control the composition of the solvent that is mixed with the steam at operation 355. At operation 365, multiphysics control system 55 uses the dielectric permittivity (∈) and the conductivity (σ) of underground deposit 20 to decide whether the total dissolved solvents (TDS) of the solvent mix have the desired composition. The dielectric permittivity (∈) and the conductivity (σ) values may be obtained by correlating the sensed temperature data (t) and sensed pressure data (P) with, for example, water phase data. The water phase data may be obtained using, for example, water phase data derived from a Mollier water phase chart.

If the solvent TDS composition is correct, multiphysics control system 55 provides the solvent at operation 370 for mixing at operation 355. However, if the TDS does not have the correct composition, the multiphysics control system 55 adjusts the composition at operation 375 before execution of operation 370.

FIG. 4 is a Mollier chart 400 showing temperature on the x-axis and pressure on the y-axis for different water phases. As shown, water is in a vapor (steam) phase in the region 403 of Mollier chart 400. Region 403 is defined by phase boundaries 405, 410, and 413. Water is in a liquid phase in region 415 of Mollier chart 400. Region 415 is defined by phase boundaries 405, 420, 425, and 430.

As shown in FIG. 3, the sensed temperature and pressure data provided by the in-situ sensor network 310 is used to determine whether the underground deposit 20 and/or well 340 is operating in region 415 of Mollier chart 400. The multiphysics control system 55 controls the electromagnetic heating system 45 to drive the well 340 and/or underground deposit 20 to an optimized portion of region 415. However, if the sensed temperature and pressure data indicate that the well 340 and/or underground deposit 20 is operating in region 403, multiphysics control system 55 executes operation 330 so that kill switch 320 effectively open circuits electromagnetic heating system 45 thereby preventing provision of further electromagnetic energy to the underground deposit 20. The multiphysics control system 55 may also open circuit kill switch 320 at operation 380 when the data provided by in-situ sensor network 310 indicates that the well 340 and/or underground deposit 20 is operating at a super-critical temperature. Although FIG. 3 shows the kill switch 320 used in connection with the application and removal of RF power to underground deposit 20, the system 15 may include other kill switches that, for example, prevent injection of the steam and solvent mix dependent on whether or not predetermined conditions are met. Further, the system 15 may use other manners of controllably applying and removing RF power to the underground deposit 20.

FIG. 5 is also a Mollier chart 500 showing temperature on the x-axis and pressure on the y-axis for different water phases. In this instance, however, the temperature data and pressure data of the Mollier chart 500 is correlated with the permittivity (∈) and conductivity (σ) of underground deposit 20. As shown, the permittivity (∈) of underground deposit 20 when system 15 is operated in the central portion of region 415 is approximately 81 F/m as determined from the Mollier data. This value varies depending on the portion of region 415 in which the system 15 is operating. The conductivity (σ) of underground deposit 20 likewise depends on the portion of region 415 in which the system 15 is operating. As shown, the conductivity (σ) of underground deposit 20 when system 15 is operated in region 415 corresponds to an induction log baseline that is determined from the Mollier data.

The well behavior model may be updated to more accurately reflect the characteristics of underground deposit 20 over time. Updates may be based on resource samples recovered from underground deposit 20 and decomposed at operation 377. Operation 377 may be executed as an online or off-line process. The decomposed materials may be separately analyzed at operation 380 and categorized, for example, as rock, water, oil, and gas. The weight percentage w % and/or the bulk volume percentage s % of each material is used at operation 382 to update the well behavior model. Additionally, or in the alternative, the sensed temperature (t) and the sensed pressure (p) from the in-situ sensor network 310 may be used at operation 382 to update the well behavior model. The updated well behavior model is employed at operation 305 to update the control parameters used to operate system 15. Still further, the updated well behavior model may be used at operation 385 to generate a substantially real-time display of the characteristics of underground deposit 20.

FIG. 6 shows one embodiment of the kill switch 320. In this embodiment, an electromagnetic signal is provided from an RF source 605 to an input of an RF driver 610. The output of RF driver 610 is provided to pressure 150 which, in turn, provides electromagnetic energy to underground deposit 20. The RF driver 610 receives its power through a power switch 615. The power switch 615 is either activated or deactivated based on the level of a switch control signal 620, which is provided at the output of a comparator 625. A reference signal 630 is provided to a first input of comparator 625 and a kill condition signal 635 is provided to a second input of the comparator 625.

When multiphysics control system 55 detects a condition in which the supply of RF power to underground deposit 20 is to cease, it directs the kill condition signal 635 to a predetermined signal level. The comparator 625 compares the reference signal 630 with the kill condition signal 635 to direct power switch 615 to an active level. When this occurs, power switch 615 drives RF driver 610 to a state in which power is removed from RF driver 610. In this state, the electromagnetic energy from RF source 605 is no longer provided to antenna 50.

FIG. 7 illustrates the impedance encountered at various portions of system 15. As shown, the RF signal generator 185 has an impedance Z_(sig), the antenna 50 has an impedance Z_(app), and the underground deposit 20 has an impedance Z_(dep). The impedance matching network 190 has an impedance Z_(match) and is disposed between RF signal generator 185 and antenna 50 with the objective of matching the impedances throughout system 15 to thereby optimize power transfer between RF signal generator 185 and underground deposit 20.

One embodiment of impedance matching network 190 includes a first input 705 and a second input 710 that are each connected to receive electromagnetic energy from RF signal generator 185. The impedance matching network 190 also may include a first output 712 and a second output 713 that are each connected to provide electromagnetic energy to antenna 50 which, in turn, provides electromagnetic energy to underground deposit 20. The first input 705 is connected to a first terminal of a first variable capacitor 715 that extends between first input 705 and node 720. A second variable capacitor 725 extends between node 720 and first output 712. A variable conductor 730 extends between node 720 and a further node 735. The node 735 is common to both the second input 710 and second output 713. The values of the first variable capacitor 715, the second variable capacitor 725, and the variable conductor 730 may be controlled by multiphysics control system 55 to alter the value of Z_(match) and achieve the desired power transfer characteristics.

As shown at operation 385 of FIG. 3, multiphysics control system 55 may display characteristics of the underground deposit 20. FIG. 8 shows one embodiment of a display system 800 that may be used to display such characteristics. The display system 800 includes multiphysics control system 55 which is connected to receive the sensed/calculated data 60 and to provide control output signals 805.

As noted in FIG. 3, multiphysics control system 55 uses and maintains a well behavior model 810 having data corresponding to characteristics of underground deposit 20. Such characteristics may include, for example, electromagnetic characteristics, fluid dynamic characteristics, resource component saturation characteristics, and the like. The multiphysics control system 55 extracts characteristic data from the well behavior model 810 to generate voxel data 815. The voxel data 815 includes volume elements representing characteristics of underground deposit 20 in three-dimensional space. Voxels for different characteristics of the underground deposit 20 overlap one another in the same volumetric space thereby allowing the multiphysics control system 55 to generate spatial maps of the characteristics of the underground deposit 20.

To reduce the amount of voxel data 815 used to generate a spatial map, the multiphysics control system 55 may provide the voxel data 815 to a blob detector 820. Various blob detection techniques may be used by the blob detector 820 including, for example, Laplacian of Gaussian blob detection, Difference of Gaussians blob detection, and Determinant of Hessian blob detection. Other blob detection techniques may also be employed, the foregoing merely being examples.

Image data generated by the blob detector 820 may be provided to an image/video processor 825. The image/video processor 825 may execute image processing operations on the data including, for example, edge detection, color mapping, pixel mapping, or the like. The processed image data may be communicated to a terminal 830 along bus 835 for display to a user. The user may provide commands through terminal 830 to select which of the characteristics of the underground deposit 20 are to be displayed on the terminal 830. Additionally, or in the alternative, the output of blob detector 820 and/or the processed image data may be provided to a further controller that implements supplementary operational features of the system 15, such as image recognition to identify portions of interest in the image. These portions of interest may be used, for example, to identify critical areas of the underground deposit 20 needing attention.

FIG. 9 is a voxel manipulation system 900 executed by the multiphysics control system 55 to generate the voxel data 815. Many of the operations executed by the voxel manipulation system 900 are used to simulate multiple physical models or multiple simultaneous physical phenomena occurring in underground deposit 20. The multiple physical models and/or multiple simultaneous physical phenomena typically involve coupled systems in which simulations for each system impacts the simulation of the other(s). The multiphysics tools used in the multiphysics control system 55 may solve partial differential equations for the coupled system to generate data corresponding to the characteristics of underground deposit 20.

The voxel manipulation system 900 of FIG. 9 includes two multiphysics processing tools that execute operations associated with two coupled systems of the underground deposit 20. In this embodiment, the voxel manipulation system 900 includes an electromagnetic tool 905 simulating the electromagnetic behavior of the underground deposit 20, and a fluid dynamics tool 910 simulating the fluid dynamic behavior of the underground deposit 20. Various partial differential equation solvers may be used to implement the electromagnetic tool 905 and fluid dynamics tool 910. For example, suitable multiphysics solvers are available from ANSYS, Inc., CPMSOL, Inc., and other providers.

The fluid dynamics tool 910 of FIG. 9 receives data from well behavior model 810. For example, well behavior model 810 may provide fluid dynamics tool 910 with pressure, porosity, saturation, viscosity, density, and mole fraction information. As discussed above, the data of well behavior model 810 may be based on an initial well behavior model as well as any updates made to that model. Updates may be based on real-time sensor data, recovered resource analysis, and/or data generated by the voxel manipulation system 900.

The fluid dynamic system modeled by the fluid dynamics tool 910 is coupled to the electromagnetic system modeled by electromagnetic tool 905. Accordingly, the fluid dynamic data provided by fluid dynamics tool 910 is dependent on values of the electromagnetic parameters modeled by electromagnetic tool 905. Using these inputs, fluid dynamics tool 910 generates specific heat voxels and mass characteristic voxels 915 modeling the underground deposit 20. The fluid dynamics tool 910 may also generate porosity voxels, saturation voxels, and density voxels 920. The density voxels 920, in turn, may be used to generate water weight percentage voxels 925.

The water weight percentage voxels 925 may be used to generate dielectric and conductivity voxels, shown collectively at 930. The conductivity voxels 930 are provided to the electromagnetic tool 905 to generate power density voxels 935. The power density voxels 935 and the mass characteristic voxels 915, in turn, are used to generate an input voxel to the fluid dynamics tool 910. Since the voxels modeling the electromagnetic system depend on voxels modeling the fluid dynamics system, and vice versa, voxel manipulation system 900 generates models for both of the coupled systems.

As shown in FIG. 9, there are several equations used to generate voxels of one characteristic and/or attribute from voxels of another characteristic and/or attribute. The voxel manipulation system 900 uses four such equations. Two equations are used at operation 940 to generate conductivity voxels 930 from water weight percentage voxels 925. In this example, conductivity voxels 930 are generated from the water weight percentage voxels 925 using Equation 1 below:

σ=5·10⁻⁴ ·w ²  (1)

where w is the water weight percentage represented by the water weight percentage voxels 925 and corresponds to a water weight percentage between 0% and 100%.

Dielectric voxels 930 may be generated from the water weight percentage voxels 925 using Equation 2 below:

∈_(Ri)(w)=a _(i) +b _(i) w+c _(i) w ²  (2)

where w is the water weight percentage represented by the water weight percentage voxels 925 and corresponds to a water weight percentage between 0% and 100%, and where a_(i), b_(i), and c_(i) are coefficients that depend upon the frequency of operation of RF signal generator 185 and the materials in underground deposit 20. FIG. 10 is a table of the coefficients at various operating frequencies for Athabasca oil sands.

The water weight percentage voxels 925 may be generated from density voxels 920 at operation 945 using Equation 3 below:

$\begin{matrix} {{w\mspace{14mu} \%} = {100*\frac{\varphi*\rho_{w}*S_{w}}{{\varphi*\left( {{\rho_{w}*S_{w}} + {\rho_{o}*S_{o}} + {\rho_{g}*S_{g}}} \right)} + {\left( {1 - \varphi} \right)*\rho_{rock}}}}} & (3) \end{matrix}$

where φ is the porosity of the oil sands, ρ_(w) is the water density, S_(w) is the water saturation, ρ_(o) is the oil density, S_(o) is the oil saturation, ρ_(g) is the gas density, S_(g) is the gas saturation, and ρ_(rock) is the rock density.

The fluid dynamics tool 910 may generate the density voxels 920 and mass characteristic voxels 915 at operation 950. At operation 950, the multiphysics control system 55 executes Equation 4 below, using the mass characteristic voxels 915 and the power density voxels 935 as inputs. The output of Equation 4 is provided as an input to fluid dynamics tool 910.

Q=c*m*ΔT  (4)

where Q is the heat energy of the voxel is, c is the specific heat of the voxel, m is the mass of the voxel, and ΔT is the temperature change of the voxel over time. Although Equation 4 does not expressly use the power density voxels 935, the Q of each voxel corresponds to the power density represented by each of the power density voxels 935.

A user may select one or more of the underground deposit 20 characteristics shown in FIG. 9 for display on terminal 830. The multiphysics control system 55 may pass the voxels of the selected characteristic to blob detector 820 to ultimately generate a map of the characteristic on terminal 830. When well behavior model 810 is frequently updated, the map on terminal 830 represents a generally real-time display of the selected characteristic.

FIGS. 11 and 12 are examples of maps that may be displayed on terminal 830. FIG. 11 is a temperature map of the underground deposit 20 at a given time, t. Region 1100 is proximate the pressure 150 and delivery pipe 35 and, as such, has the highest temperature. The temperature decreases with distance from region 1100. FIG. 12 is an oil saturation map of underground deposit 20 at a given time, t. Again, region 1100 is proximate the pressure 150 and delivery pipe 35 and, as such, has the highest oil saturation. As noted above, other characteristics of underground deposit 20 may also be presented to a user on such maps.

Although preferred embodiments of the invention have been described using specific terms, devices, and methods, such description is for illustrative purposes only. The words used are words of description rather than of limitation. It is to be understood that changes and variations may be made by those of ordinary skill in the art without departing from the spirit or the scope of the present invention, which is set forth in the following claims. In addition, it should be understood that aspects of the various embodiments may be interchanged either in whole or in part. Therefore, the spirit and scope of the appended claims should not be limited to the description of the preferred versions contained herein. 

1. A control system for use in extracting hydrocarbons from an underground deposit comprising: an electromagnetic heating system configured to heat the underground deposit to facilitate fluid flow of a resource through the underground deposit; and a processor configured to control the electromagnetic heating system in response to temperature data and pressure data for the underground deposit, wherein the processor is further configured to correlate the temperature data and the pressure data with predetermined water phase characteristics to control the electromagnetic heating system so that substantially all water in the underground deposit is maintained in a liquid state.
 2. The control system of claim 1, wherein the processor is configured to use correlate the temperature data and the pressure data with Mollier water phase temperature and pressure correlations to determine water phase characteristics within the underground deposit.
 3. The control system of claim 1, wherein the electromagnetic heating system comprises: an electromagnetic energy source having variable frequency characteristics and variable power output characteristics; and an antenna configured to provide electromagnetic energy from the electromagnetic energy source to the underground deposit.
 4. The control system of claim 3, wherein the processor is further configured to control the variable frequency characteristics and the variable power output characteristics of the electromagnetic energy source in response to electromagnetic parameters of the underground deposit.
 5. The control system of claim 4, wherein the electromagnetic parameters comprise permittivity and conductivity of the underground deposit.
 6. The control system of claim 5, wherein the processor is configured to determine the permittivity and the conductivity using the temperature data and the pressure data.
 7. The control system of claim 1, further comprising a display system configured to provide an image of at least one attribute of the underground deposit during extraction of the hydrocarbons.
 8. A control system for use in extracting hydrocarbons from an underground deposit comprising: an electromagnetic energy source; an antenna configured in the underground deposit to heat the underground deposit in response to electromagnetic energy received from the electromagnetic energy source; an impedance matching network configured to communicate electromagnetic energy between the electromagnetic energy source and the antenna; a temperature sensor configured to provide a sensed temperature of the underground deposit; a pressure sensor configured to provide a sensed pressure of the underground deposit; and a processor configured to control electrical parameters of the electromagnetic energy source and the impedance matching network in response to permittivity and conductivity of the underground deposit, wherein the processor is further configured to correlate the sensed temperature and the sensed pressure with one another to determine the permittivity and the conductivity of the underground deposit.
 9. The control system of claim 8, wherein the processor is configured to correlate the sensed temperature and the sensed pressure with predetermined water phase characteristics to control the electromagnetic energy source so that substantially all water in the underground deposit is maintained in a liquid state.
 10. The control system of claim 8, further comprising a display system configured to provide an image of at least one attribute of the underground deposit during extraction of the hydrocarbons.
 11. A system for providing an image of an underground deposit during extraction of hydrocarbons comprising: a multiphysics control system having a well behavior model including data corresponding to characteristics of the underground deposit, a multiphysics processing tool configured to generate characteristic data for the underground deposit using the well behavior model and a model of a coupled system of the underground deposit, a voxel generator configured to generate voxels for the underground deposit using the characteristic data generated by the multiphysics processing tool; a display; and an image processor configured to provide image data to the display, wherein the image data is derived from the voxels generated by the voxel generator.
 12. The system of claim 11, wherein the data of the well behavior model comprises electromagnetic characteristic data.
 13. The system of claim 11, wherein the data of the well behavior model comprises fluid dynamic characteristic data.
 14. The system of claim 11, wherein the data of the well behavior model comprises resource component saturation characteristic data.
 15. The system of claim 11, further comprising a blob detector configured to execute blob detection operations on the voxels generated by the voxel generator.
 16. The system of claim 15, wherein the blob detection operations are selected from a group of blob detection operations consisting of Laplacian of Gaussian blob detection, Difference of Gaussians blob detection, and Determinant of Hessian blob detection.
 17. The system of claim 11, wherein the well behavior model is configured to be updated using real-time temperature and real-time pressure data.
 18. The system of claim 11, wherein the well behavior model is configured to be updated using underground deposit data derived from an extracted resource analysis.
 19. The system of claim 11, wherein the image processor is configured to generate a color map of one or more of the characteristics of the underground deposit.
 20. The system of claim 19, wherein the color map corresponds to a state of the one or more of the characteristics of the underground deposit occurring at a given time.
 21. The system of claim 11, multiphysics processing tool correlates fluid dynamic characteristics of the underground deposit with electromagnetic characteristics of the underground deposit to model the coupled system.
 22. The system of claim 21, wherein the multiphysics processing tool comprises a fluid dynamics tool and an electromagnetic tool.
 23. The system of claim 22, wherein the fluid dynamics tool is configured to generate specific heat voxels and mass characteristic voxels for the underground deposit.
 24. The system of claim 22, wherein the fluid dynamics tool is further configured to generate porosity voxels, saturation and voxels, and density voxels.
 25. The system of claim 11, wherein the voxel generator is configured to generate conductivity voxels from water weight percentage voxels.
 26. The system of claim 11, wherein the voxel generator is configured to generate dielectric voxels from water weight percentage voxels.
 27. The system of claim 11, wherein the voxel generator is configured to generate water weight percentage voxels from density voxels. 